Oxygen vacancies at titanate interfaces: 
two-dimensional magnetism and orbital reconstruction 



N. Pavlenko,^'2'El T. Kopp^ E.Y. Tsymbal,^ J. Mannhart,"^ and G.A. Sawatzky^ 

^ Center for Electronic Correlations and Magnetism, Theoretical Physics III, 
Institute of Physics, University of Augsburg, 86135 Augsburg, Germany 
'^Center for Electronic Correlations and Magnetism, Experimental Physics VI, 
Institute of Physics, University of Augsburg, 86135 Augsburg, Cermany 
^Department of Physics and Astronomy, Nebraska Center for Materials and Nanoscience, 
University of Nebraska, Lincoln, Nebraska 68588-0299, USA 
^Max Planck Institute for Solid State Research, 70569 Stuttgart, Cermany 
^Department of Physics and Astronomy, University of British Columbia, Vancouver, Canada V6T1Z1 

We show that oxygen vacancies at titanate interfaces induce a complex multiorbital reconstruction 
which involves a lowering of the local symmetry and an inversion of t2g and Cg orbitals resulting in 
the occupation of the Cg orbitals of Ti atoms neighboring the O vacancy. The orbital reconstruction 
depends strongly on the clustering of O vacancies and can be accompanied by a magnetic splitting 
between the local Cg orbitals with lobes directed towards the vacancy and interface dxy orbitals. The 
reconstruction generates a two-dimensional interface magnetic state not observed in bulk SrTiOs. 
Using generalized gradient approximation (LSDA) with intra-atomic Coulomb repulsion (GGA+U), 
we find that this magnetic state is common for titanate surfaces and interfaces. 

PACS numbers: 74.81.-g,74.78.-w,73.20.-r,73.20.Mf 



I. INTRODUCTION 

Charged impurities can change the electronic prop- 
erties of insulating materials qualitatively. For exam- 
ple, doping of SrTiOa with Nb leads to semiconducting 
properties and eventually to a superconducting transi- 
tion below 0.3 K^. The recent discovery of the metal- 
lic, superconducting and magnetic states at the in- 
terface between the bulk insulators SrTiOs (STO) and 
LaA103(LA0) triggered an intense exploration of elec- 
tronic reconstruction and the role of impurities in the 
formation of a conducting and sometimes magnetic state 
at these interfaces^"—. 

Two distinct doping mechanisms, controlled by the 
growth conditions of the samples, appear to be responsi- 
ble for the formation of interface electron carrier^—. 
The intrinsic polar charge doping in the heterostruc- 
tures prepared at high oxygen pressures (10~^ mbar) 
suggests doping with an upper carrier density limit of 
3 • 10^^ cm~^ which corresponds to 0.5 electrons per in- 
terface unit cell (uc) and compensates the interface polar 
discontinuity. The resulting two-dimensional electron liq- 
uid forms in a nm-thick interface layer— i^*^ — and can 
be a possible source for superconducting properties tuned 
by gate electrostatic fieid^^. The extrinsic mecha- 
nism of charge doping due to oxygen vacancies domi- 
nates in the LAO/STO samples grown at lower oxygen 
pressures (10~^ mbar) reported in Refs. Q and 0. Re- 
markably, the oxygen reduced heterostructures are char- 
acterized by higher charge densities of 10^^ — 10^^ cm~^ 
with the charge spatially extended into the /im-thick in- 
terfacial layer inside the STO^. Recent cathode and 
photoluminescence experiments^ provide direct evidence 
for oxygen vacancies in the STO-substrate. The oxy- 
gen vacancies can diffuse in STO with activation energy 



of 0.75 eV, which is comparable to the activation energy 
1 eV of the vacancy migration through the LAO layers^. 
At each oxygen vacancy in the Ti02 layer, two electrons 
are weakly bound to retain charge neutrality of a config- 
uration in which the two titanium ions in the dimer share 
two electrons. Due to these moderately bound electrons 
the oxygen vacancies act as n-type dopants. We note 
though that in STO, these O- vacancy induced, weakly 
bound electrons do not contribute to resolve the polar 
discontinuity problem since they simply are there to pro- 
duce charge neutrality upon removing the oxygens. The 
0-vacancy induced electrons are bound to the effective 
2+ charge of the vacancy. Below we argue that these elec- 
trons enhance the interface conductivity considerably. 

A critical issue related to the electron concentration is 
the electric conductivity of LAO/STO heterointer faces. 
The experimental measurements show that the conduc- 
tivity of n-type LAO/STO bilayers increases with low- 
ering the oxygen partial pressure po^ and exhibits an 
abrupt jump from 10~^ to 10^ /n when pos changes 
from 10~^ to 10~^ mbar— i^. The high conductivity shown 
by samples grown at low oxygen pressures is caused by a 
large amount of oxygen vacancies which act as electron- 
donors in such samples^' 

The local electronic state of the vacancy-containing 
transition metal oxides can exhibit extraordinary prop- 
erties not observed in chemically stoichiometric materi- 
als. Recent first principle studies demonstrate a vacancy- 
related magnetic exchange splitting in nonmagnetic ma- 
terials like CaO or LAO/STO^^^^^ For LAO/STO, cal- 
culations of the vacancy formation energiesi^i^ suggest a 
predominant positioning of 0-vacancies in the top AIO2 
surfaces. The energy barrier for transport of 0-vacancy 
from the interface Ti02 to the surface AIO2 layer strongly 
depends on the vacancy concentratio n^^i^^ . 
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Recent STM, cathode luminescence studies and con- 
ductivity measurements provide strong support for clus- 
tering of the oxygen vacancies in STO^i^. Moreover, the 
STM experiments and DFT-\-U studies suggest the possi- 
bility that oxygen vacancies form pieces of linear chains, 
which appear to be more energetically favourable^^ as 
compared to the isolated vacancies^ —. 

In this work, we explore the electronic structure of the 
first titanate surface/interface layers considering different 
configurations and concentrations of oxygen vacancies. 
In bulk SrTiOs, the octahedral crystal ligand field, co- 
valency and Coulomb repulsion result in higher energies 
of Ti eg states as compared to t2g states with a split- 
ting of about 2.5 eV— — . Here we find that in contrast 
to the bulk stoichiometric SrTiOs, each Ti-O^-Ti dimer 
contains missing hybridization links which strongly lower 
the energy of the d orbitals with lobes pointing to the va- 
cancy. As a result, we obtain a splitting of the doubly de- 
generate Cg orbitals and the triply degenerate t2g orbitals 
into orbital doublets and singlets. In the pure material, 
these e^-like antibonding orbitals are strongly pushed 
up due to the strong hybrization with the lower-energy 
O 2p orbitals. In the oxygen-reduced LAO/STO, the 
missing local covalent bonding and the local symmetry 
lowering results in a new kind of orbital reconstruction 
at titanate interfaces: a shift and partial occupation of 
the vacancy-directed Ti e^-orbitals accompanied by their 
magnetic splitting and mixing with the 3dxy states. Our 
calculations of SrTiOa-surfaces and LAO/STO interfaces 
demonstrate the universal character of the orbital recon- 
struction of the titanates due to surface/interface oxygen 
vacancies. In the linear clusters of vacancy stripes, in 
which each Ti is accompanied by two 0-vacancies the Cg 
orbital state pointing along the Ti vacancy direction, is 
further strongly lowered. In these Cg states, both spin 
states are occupied, although a small amount of spin po- 
larization is also present. The high spin polarization oc- 
curs predominantly in the t2^-states. 

We also show that the electronic states of the oxygen- 
reduced STO and LAO layers are physically different. 
The vacancies in Ti02 layers release extra 3d electrons 
which are confined to the quasi-two-dimensional layer of 
the electron liquid at the LAO/STO interface. In con- 
trast, the oxygen vacancies in the LAO overlayer pro- 
duce local electronic states of a mixed sp character which 
are redistributed between the surface AIO2 and interface 
Ti02 layers resulting in the confinement of electrons in 
the parallel, surface and interface 2D layers with different 
effective masses, a concept which has been recently dis- 
cussed in Refs. |4l| and'42. We find that for small concen- 
trations of oxygen vacancies cy < 1/8 at the AIO2 sur- 
face, the vacancy released electrons are transferred to the 
LAO/STO interface and compensate the interface polar- 
discontinuity, which also leads to an insulating surface 
state. 



II. OXYGEN VACANCIES AT LAO/STO 
INTERFACES 

We use supercells in density functional theory to ex- 
plore the charge and spin state around oxygen vacancies 
at LAO/STO heterointer faces. The supercells contain a 
4 uc thick LAO overlayer deposited on a STO layer of a 
thickness between 1 uc and 6 uc. The LAO-STO-LAO 
structures are separated by a 13 A thick vacuum sheet. 
In the supercell, oxygen vacancies can be located in the 
interface Ti02 layer or in one of the A102-planes of LAO. 
Cells with three types of oxygen vacancy configurations 
in MO2 (M=Ti, Al) are sketched in Fig.I^a), (b) and (c). 
Configuration (a) contains oxygen vacancies in Ti-O^-Ti 
dimers, with each vacancy located between two Ti atoms 
in a (\/2 x ^/2) uc so that each Ti in the Ti02 has exactly 
one nearest vacancy. We refer to this as a Ti dimer struc- 
ture. In contrast, the configuration (b) is introduced by 
removing the oxygen atom 0(a/2, b/2) in the center of 
the (2x1) M2 04-plaquette, a configuration which results 
in stripes of vacancies in the ^/-direction. The doping level 
of 0-vacancies in configurations (a) and (b) corresponds 
to cv = 1/4 (25%) of vacancies per each four oxygen 
atomic sites. The lowest density of 0-vacancies (one va- 
cancy per 8 oxygen sites) is represented by the config- 
uration (c) in Fig. [TJ which as in configuration (a) has 
Ti-O^-Ti dimers but the density of vacancies is a factor 
of two smaller. 

The density functional calculations were performed 
using the Generalized Gradient Approximation (GGA) 
in the Perdew-Burke-Ernzerhof pseudopotential 
implementation^ in the Quantum Espresso (QE) 
package^^. We use a kinetic energy cutoff of 640 eV, and 
the Brillouin zone of the 106- to 166-atom supercells 
is sampled with 5x5x1 to 9x9x1 /c-point grids. An 
increase of the k-point mesh from (5x5x1) to (7x7x1) 
leads to a negligibly small change of the total energy 
by 0.005 Ry and to an increase of the Ti magnetic 
moments by small values of about 0.05 /i^ for the case 
that 0-vacancies are present. In our calculations we 
account for a local Coulomb repulsion of Ti 3d electrons 
by employing a GGA+/7 approach with Uti = 2 eV— . 
The supercells have been structurally relaxed by a 
combination of the optimization procedures of the full 
potential WIEN2k-package and the pseudopotential QE 
package^i^. The in-plane lattice constants have been 
fixed to their room-temperature bulk- STO cubic values 
(<^STO = ^STO = 3.905 A). Although low-temperature 
SrTiOs is tetragonal with the lattice parameters 
a = 3.896A and c = 3.899 A, in our studies we assume 
that the small difference between the lattice constants 
should not seriously affect the state of oxygen vacancies. 



A. Vacancy stripe configurations 

We start with the confiuration with the largest degree 
of clustering of 0-vacancies in our models which is the 
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FIG. 1. Schematic view of the SrTiOs/LaAlOs heterostmc- 
ture. The supercell contains a 4 uc thick LaAlOs layer de- 
posited on a 2.5 uc thick SrTiOs slab. The full supercell 
consists of two symmetric parts of the depicted structure 
and a vacuum layer of 13 A. The structures on the right 
side show Mn02n (M=Ti, Al)-plaquettes (a) with one elim- 
inated 0(2a,a) atom (dimerized vacancy), (b) with elimi- 
nated O(0.5a,0.5a) (chained vacancy) and (c) with eliminated 
O(0.5a, 0.25a), generated for the study of the systems with 
O- vacancies. The position of an O- vacancy is identified by a 
red dashed circle. 



vacancy stripes of type (b). Fig. [2] presents the spin po- 
larized projected densities of 3d states for the Ti(0.5a,0) 
atom nearest to the 0-vacancies of type (b) and for the 
Ti(0,0) at the corners of the Ti2 04 plaquette in Fig.lTJb) 
located still in the interface layer but furthest away from 
the O vacancy stripe. For the Ti(0.5a,0) atom (Fig.[2fa)), 
we find a strong "splitting" of the 3(^3^2 and 3(1^2 
states with the 3d^y2_^2 being almost completely occu- 
pied with close to 2 electrons and the 3dz2-x'^ almost 
completely empty and at very high energies above Ef. 
The extremely narrow eg DOS peaks are located in the 
energy window {Ep — 2.3 eV^Ep + 0.5 eV) relative to the 
Fermi level. This is in contrast to the DOS of the Ti(0,0) 
which is more distant to the oxygen vacancy (Fig. [2fc)). 
At Ti(0,0) the empty eg states are located about 3.5 eV 
above Ef and 2 eV above the t2g states, the latter being 
partially occupied by the electrons generated due to the 
polar interface. 

To better understand the nature of the sharp eg peaks, 
we plot in Fig.[2](b) the 3(i3^2_^2 density of states and see 
that actually the DOS in the sharp peaks is almost solely 
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FIG. 2. Spin-polarized orbital-projected densities of states for 
the chained vacancy configuration of Fig. [TJb): (a) DOS for 
the interface Ti(0.5a,0) next to the O-vacancy(0.5a,0.5a) in 
the Ti02 layer; (b) the Ti(0.5a,0) densities of states projected 
on the rotated 3d orbitals. Here the reconstructed Cg orbital 
state corresponds to 3d^^2_r2 along the Ti-O^-Ti and the 3dxz 
corresponds to the old 3dxy state, (c) The densities of states 
for the interface Ti(0,0) more distant from the O-vacancy. 
The considered supercells comprise 4 uc thick LaAlOs layers 
and a 4 uc thick SrTiOs layer. 



in an eg 3d^y2_^2 orbital with the lobes pointing along 
the Ti-Ov direction, which implies a strong confinement 
of the vacancy-released electrons between the nearest Ti 



atoms. This strong preference of the 3d' 



orbital oc- 



cupation is clearly seen also in the spatial charge density 
distribution in FigSJ^b). The very sharp peaks in the eg 
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density of states is a result of the quasi-one-dimensional 
character of the band structure of these mainly 8(^3^2 
composed bands. For charge neutrality we need to acco- 
modate two electrons per Ti in the dimerized chain which 
aside from the polar induced charge indicates that the Ti 
ions in the chains are formally divalent with nearly two 
electrons each. We see indeed from Fig. EJb) that both 
the spin up and down 8(^3^2 states are nearly fully 
occupied. These states are therefore rather magnetically 
inert. The empty 8(^3^2 _^2 stripe projected DOS, also 
clearly visible in Fig.[2jb), results from the other Ti ions 
which do not neighbor oxygen vacancies. 

The extreme splitting of the Cg states 3dsy2_^2 di- 
rected along the Ti- vacancy direction reflects the inter- 
face orbital reconstruction, a phenomenon recently ob- 
served at the interfaces between (Y,Ca)Ba2Cu307 and 
Lao.67Cao.33Mn03^^ and predicted for interfaces between 
YBa2Cu306 and STO^^. For the oxygen deflcient 
LAO/STO interfaces, the interface orbital reconstruction 
results in the change of the chemical valence state of the 
Ti atoms close to the O- vacancies. For high enough va- 
cancy concentrations this strong valence change might 
be visible in experiments like those of Ref. |8|. More- 
over, the conclusion of Ref. [H about the in-plane mixed 
3dxy/Sdx2_y2 character of the interface Ti electrons 
strongly supports the notion of the splitting of Cg states 
and partial occupation of the in-plane vacancy-directed 
Cg orbital found in our calculations. 
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FIG. 3. Orbital-projected charge density profiles along the 
[001] direction. Here the 0-vacancy of type (b) is located 
in the TiO 2 interfacial layer in the supercell containing 4 uc 
thick LaAlOs layers and a 4 uc thick SrTiOs layer. Note the 
large scale change between the upper and lower panels. 

To estimate the chemical valence of the Ti atoms, we 
have calculated the 2D electron charge density in two dif- 
ferent energy windows below the Fermi level correspond- 
ing to occupied Cg and t2g orbitals, and performed the in- 
tegration in {x^y) planes across the interface. The results 



are presented in Fig. [H The integration of the Cg charge 
proflle along the z-axis gives a total of about 1.9 elec- 
trons, which implies that the excess vacancy-produced 
charge is almost fully transferred to the Cg orbitals and 
induces a change of the chemical valence of the near- 
est Ti atoms by —1.9. In contrast to the stoichiometric 
STO, the local e^-states originating from the vacancies in 
the interface TiO 2 layers of the STO are placed well be- 
low the Fermi level. The clustering of the vacancies into 
stripes at the Ti02/LaO interface leads to a strong con- 
centration of the vacancy-released electrons due to the 
shift of the fllled Cg subbands by more than 1 eV be- 
low the Fermi level. Integration of the conducting {t2g) 
charge proflle shown in Fig. [3] reveals that the charge of 
the 0.8 t2g electrons in the Ti02 layers near the interface 
is compensated by the charge of the holes present at the 
AIO2 surface of the LAO fllm. 

In contrast to the strong spin polarization of the t2g 
states observed in Fig. [21 the occupied Cg states are al- 
most unpolarized. The integration of the sharp Cg peaks 
below the Fermi level gives a negligibly small value of the 
spin polarization of about 0.001 /i^, which implies the ex- 
istence of the orbital separation of the spin and charge de- 
grees of freedom in the vacancy stripe conflgurations: the 
vacancy-released charge is localized in the spin-nonpolar 
e^-orbitals, whereas only the t2g polarity-induced intrin- 
sic charge contributes to the magnetization. The con- 
tribution from different orbitals to the magnetization of 
various vacancy conflgurations is presented in Table H 

Fig. [4] shows the spatial distribution of the charge den- 
sity generated in the energy windows below the Fermi 
level for the fllled t2g (Fig.HJa)) and Cg states (Fig. [Kb)). 
One can see that the polarization-generated t2g elec- 
trons mostly occupy the Ti 3dxy and O 2py orbitals 
which are hybridized and do not contain any contribu- 
tion from the oxygen vacancy state. In contrast, the 
vacancy-generated electronic charge in Fig. [H^b) has a 
mixed 3dx2-y2 /3dsz^-r^ character (or 3dsy2_r^ charac- 
ter in the projection along the Ti-O^-Ti). 

Due to Coulomb attraction with the effective 2+ 
charge of the nearest oxygen vacancy, the excess elec- 
tron charge is trapped in the in neighboring Ti sites and 
because of the change in local Cg -O 2p hybridization the 
charge is mainly found in Ti orbitals of 8(^3^2 _^2 charac- 
ter. These are placed around the center of the 0-vacancy, 
similar to the bulk STO with point defect^. 



B. Ti dimer-type O- vacancies in 
(a/2 X A/2)-plaquettes 

The strong localization of the excess electrons is al- 
ready noted in the oxygen-reduced bulk STO, where the 
tendency for the formation of local Ti states is connected 
to the degree of clustering of the oxygen vacancies^. 
Similar to bulk STO, a weaker clustering in LAO/STO 
heterointerfaces reduces the energy gap between the fllled 
and empty Ti Cg states to less than 1 eV, which is demon- 
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FIG. 4. Electron density maps for vacancy stripes as shown 
in Fig. 1(b). (a) t2g and (b) Cg electron states for the inter- 
face Ti02 layer with the O(0.5a,0.5b)-vacancy in the supercell 
containing 4 uc thick LaAlOs layers and a 4 uc thick SrTiOs 
layer. The contour lines in (a) indicate the pd-hybridization 
of the polar interface-induced electrons in Ti02 layer. 



strated in Fig.[5]for the interface Ti02 with Ti dimer-type 
0-vacancy configurations (a). 

In Fig. [5^, the filled midgap Cg states in the energy 
window {Ef — 1.2 eV; Ep — O.! eV) are strongly confined 
on Ti Cg orbitals centered around the vacancies although 
the occupied Cg subband is broad as compared to the 
extremely narrow Cg DOS-peaks in the striped configu- 
rations. The relaxation of the atomic positions leads to 
an elongation of the dimer Ti-O^-Ti by 0.1 A due to the 
shifts of the Ti atoms by about 0.05 A outward from the 
0-vacancy. The resulting DOS for the relaxed structure 
is presented in Fig. [Sb. The comparison with the DOS 
for the unrelaxed structure (top panel) shows a crucial 
effect of the structural relaxation of atoms around the 
0-vacancy. The relaxation leads to an additional mag- 
netic splitting of the Ti Cg states and to a decrease of the 
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FIG. 5. Spin-polarized orbital-projected densities of states 
for the interface Ti in the interface Ti02 layer with dimerized 
configurations (a) of oxygen vacancies. Here the top panel (a) 
represents the DOS of a structure with unrelaxed Ti atoms 
fixed in the center and in the corners of the Ti204 plaque- 
tte; the bottom panel (b) presents the 3d states of the re- 
laxed structure with the Ti atoms shifted by 0.05 A outward 
from the 0-vacancies which corresponds to the elongation of 
a vacancy- containing Ti-0^;-Ti dimer by 0.1 A. 



gap between the electronically occupied and the empty 
spin-up (majority) Cg states. The occupied Cg states are 
located in the energy window {Ep — 1.5 eV; Ep) and 
hybridize with the t2g states. As a consequence, the elec- 
tronic state of the Ti-O^-Ti dimer is half-metallic with 
conducting electrons of mixed eg-t2g character released in 
part by the 0-vacancy and from the polar induced elec- 
tronic reconstruction involving mainly t2g states. The 
mixed eg-d^y character of the electrons in the Ti-O^-Ti 
dimer reflects itself in the interface charge-density plot in 
Fig. [6l Due to the additional magnetic splitting of the eg 
states, the magnetic moment of Ti in the relaxed struc- 
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FIG. 6. Charge density plot for an interface Ti02 layer with 
one 0-vacancy at (0.75a, 0.75a) (a = 5.523 A) for each four O 
in the Ti204 plaquette in the supercell that containins two 4 
uc thick LaAlOs layers and one 4 uc thick SrTiOs layer. The 
plot has been obtained by calculating the electron densities 
in the energy window {Ef — 1.5 eV; Ef) which corresponds 
to the mixed eg-\-t2g states of Ti. 



ture is strongly enhanced to M^^f = 0.475 /i^, as com- 
pared to the unrelaxed structure with M^f^^^ = 0.15 i^b 
(see Table IH). 



C. Ti dimer-type O-vacancies in (2 x 2)-plaquettes 

To analyze lower concentrations of oxygen vacancies, 
we studied larger supercells which contain (2 x 2) M4O8- 
plaquettes in the (x^y) planes (Fig. [He)). The elim- 



TABLE I. Calculated magnetic moments of Ti ions near- 
est to the 0-vacancy and spin polarization character for 
(LaA103)4/(SrTi03)4 heterostructures with different config- 
urational state n — (i) of 0-vac. Here n is the index of the 
M02-layer in STO or LAO with top surface/interface layer 
corresponding to n = 0, and i = a, b,c is the type of vacancy 
configuration. 
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0-(a) (STO)-relaxed 


t2g + eg 


0.475 


0-(a) (STO)-unrelaxed 


t2g 


0.15 


0-(b) (STO) 


t2g 


0.34 


0-(c) (STO) 


t2g-\-eg 


0.47 


l-(b) (STO) 


t2g 


0.12 


2-(b) (STO) 


t2g 


0.07 


3-(b) (STO) 


t2g 


0.006 


0-(b) (LAO) 


t2g 


0.22 


0-(c) (LAO) 


t2g+eg 


0.56 



o 
o 



O 



0.5 



\ 




-0.5 - 



o 

in 
o 



total (up) 
total e (dn) 











— \ — 


— \ — 


— \ — 






H 



in 
O 



3d^^ 
3d 



in 
O 
O 



(a) 



0.5 



-4 -2 
Energy E-E (eV) 






v 

-0.5 - 



O 

O -1 



T 



T 



total (up) 
total (dn) 

H \ 



i«PS*:5g^^»-■ 



3d^, 
3d 




(b) 



-4 -2 
Energy E-E (eV) 



FIG. 7. Spin-polarized orbital-projected densities of states 
for the interface Ti in the interface Ti02 2 x 2-plaquette con- 
figuration with 1/8-concentration of oxygen vacancies of tpe 
(c). The positions of Ti, Al, La and O atoms have been fully 
relaxed in the (x, ^)-planes and in the z-direction. The top 
(a) and bottom (b) plots represent the DOS for the Ti atoms 
nearest to the 0-vacancy and the more distant corner Ti atom 
of the Ti408-plaquette. 



ination of one oxygen atom in the interface (Ti02)4- 
plaquette corresponds to cy = 1/8-concentration or 
1.5 • 10^^ cm~^ density of vacancies homogeneously dis- 
tributed in the interface Ti02 layer. The orbital pro- 
jected Ti 3d DOS for the vacancy configuration of type 
(c) is displayed in Fig. [71 As shown, even the smaller 
concentration of oxygen vacancies still causes a splitting 
of the Cg orbitals of Ti atoms in the close proximity of 
the vacancy. Similar to the other dimerized configura- 
tions, the occupied Cg and t2g states are located in the 
same energy window {Ef — 1.5 eV; Ef) which implies a 
mixed eg-t2g character of the vacancy-generated electron 
states. Fig. [8l^a) presents the conducting electron den- 
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FIG. 8. Profiles of the density of mobile charges along the 
[001] direction with Cv = 1/8-concentration of oxygen vacan- 
cies of type (c). The two profiles correspond to the 0-vacancy 
placed (a) in the Ti02 interface and (b) in the AIO2 surface 
layer of (2 x 2)-LA0/ST0 superceh. 
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FIG. 9. Profiles of x^/- averaged electrostatic potentials along 
the [001] direction with cv = 1/8-concentration of oxygen 
vacancies of type (c) at the interface. The black and red 
profiles correspond to the stoichiometric structure and the 
structure with the 0-vacancy in the Ti02 interface of (2x2)- 
LAO/STO supercell, respectively. 



sity profile across the interface calculated by integrating 
the planar charge density generated in the energy range 
{Ef — 1.5 eV; Ep). The location of the oxygen vacancies 
in STO leads to an insulating state of the A102-surface. 
In Fig. [3 the density plots are restricted to the interface 
region and are zero beyond 3A from the interface due to 
insulating character of the top AIO2 layer. Moreover, the 
integration of the charge profile in Fig. [8fa) gives exactly 
two electronic charges, in accordance with the charge re- 
leased by the oxygen vacancy, which implies the absence 
of the polarity-induced interface electrons and a suppres- 
sion of the polar character of the LAO/STO interface. 
The reason for this at first glance surprising effect is dis- 
cussed in the next section. 



D. Vacancy-enhanced L AO-critical thickness 

An examination of the atomic positions in the struc- 
tures with interface 0-vacancies reveals a considerable 
increase of the atomic distortions near the vacancy in 
both the interface Ti02 and in the LaO layers, presented 
in Table im These distortions originate from the repulsion 
between La ionic charges and the nearest positive charged 
0-vacancy. Table [III demonstrates a drastic increase of 
the La-0 and Ti-0 buckling which is related to an elon- 
gation of the local LaO bonds by about 0.2 — 0.3 A in the 
LAO layer. The distortions in the LaO-layers have an ex- 
tended character along the [001] -direction and are signif- 
icant at the interface as well as in the surface LaO-layers, 
which induces additional dipole moments antiparallel to 
the LAO layer polarization Plao = — A'LAoCLAoe/2 = 
-7.578 eA (here 7VlaoClao/2 = 7.578 A is the thick- 
ness of the LAO layer). The dipole moments contribute 
to the compensation of the polarization field and, as 
a result, prevent the electronic reconstruction on ac- 
count of a compensating polar charge. In the structure 
with Cv = 1/8 vacancy concentration in the interface 
Ti02, the additional vacancy-induced dipole polarization 
is APTi02 ^ 0-9 eA in the interface Ti02 layer and 
^PLaO = 0-86 eA in the subsurface LaO layer. This 
leads to the compensation of Plao by the sum of antipar- 
allel contributions (1) APlao (stoich) originating from 
the distortions in the stoichiometric structure and (2) 
APvac caused by the enhanced distortions in the LAO 
layers due to 0-vacancies in the Ti02 layer. The preven- 
tion of the electronic reconstruction and intrinsic doping 
due to enhanced distortions in the LAO-layer implies a 
vacancy-driven enhancement of the LAO critical thick- 
ness below which the polar problem can be solved by 
atomic displacements in the LAO layer— 

Fig. [9] presents the microscopic x?/-averaged and 
macroscopically averaged electrostatic potentials calcu- 
lated in the stoichiometric LAO/STO with 4 u.c. thick 
LAO layer, and in the structure with cy = 1/8 oxygen 



TABLE II. Atomic displacements (in A) in the interface Ti02 
and different (LaO)n (n = 1, . . . , 4) planes of STO/LAO with 
one 0-vacancy of type (c) in the interface Ti02 {cv = 1/8)- 
Here the buckling of each layer is determined by the maximal 
atomic distortions in the corresponding MO-layer (M=Ti,La) 
which are defined as A^^mo = A^m — A^o- The values Ao 
correspond to the buckling in the stoichiometric vacancy-free 
structure. 



(M02)n 


A^MO 


Ao2;mo 


A- Ao 


(Ti02)i 


0.21 


-0.1 


0.31 


(LaO)i 


0.33 


0.03 


0.3 


(LaO)2 


0.28 


0.05 


0.23 


(LaO)3 


0.28 


0.07 


0.21 


(LaO)4 


0.36 


0.08 


0.32 
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38 



FIG. 10. Schematic view of a Tii-0^-Ti2-cluster that models 
an oxygen vacancy with the corresponding reconstruction of 
the 3d electron levels. Here the bare unrenormized energies 
Sp - -10.199 eV and Sd = -6.01 eV; Agt = 2.7 eV, Ai 



3.2 eV, and A 



gt 



-2 eV. 



vacancies in the interface Ti02-layer. The macroscopic 
potential is calculated by the macroscopic averaging pro- 
cedure proposed in Ref. "51 . As compared to the stoichio- 
metric system with the macroscopic electric field Egt = 
0.17 eV/A screened by the polar charge, the macroscopic 
electric field in the structure with the interface vacancies 
is screened by the field APvac^ generated by antipolar 
distortions and approaches the value Eye = 0.13 eV/A. 
For the obtained internal field 0.13 eV/A, the critical 
LAO-thickess sufficient for the dielectric breakdown of 
LAO/STO with the STO band gap 3.2 eV^i in the struc- 
ture with the interface 0-vacancies is enhanced up to 
24.6 A which corresponds to a 7 u.c. thick LAO layer. 



III. A MODEL FOR THE Cg LEVEL SPLITTING 

The vacancy-induced Ti 3d orbital reconstruction can 
be understood from the analysis of the local bonding for 
a two-electron state of a Ti-O^-Ti-cluster (Fig. [TQ|) . 

We focus on the contribution of the covalent Ti3(i-02p 
bonding to the diagonal and exchange contributions to 
the electron energy levels which determine the magnetic 
state (paramagnetic singlet or ferromagnetic triplet) of 
the two-electron cluster. 

In the stoichiometric SrTiOs, the two-electron energy 
contains a pd-hybridization term due to the overlap be- 
tween the 3(i-states ot Ti atoms and 2j9-states of the cen- 
tral oxygen atom. In a TiOe-octahedron, such a covalent 
contribution Vpd produces shifts of the one-electron 3d 



Here St = -6.258 eV and Se = -5.638 eV are the Ti 
3d ionization energies plus Madelung potential, renor- 
malized by the electrostatic shifts due to the cubic field. 
The corresponding energies for the O 2p states are: 
£11 = -10.519 eV and ei_ = -10.039 eV. In Eq. Sa = 

x^y^z) and s'^ -- 



^(3 ^a^% 



coskaa {a = x,y,z) and = y " 
With the pd-covalency parameters Vpda ^ 2.1eV and 
VpdT^ ~ 0.84eV defined in Refs. IH and ST 



we can es- 



timate the splitting between the anti-bonding eg- and 
t2^-states in the F-point of the Brillouin zone^^i"— : 



^gt 



e{E)-s{T2) 



pda 



PCITV 



Ao 



2.7 eV, (2) 



where Vpda = V^Vpda, Vpdi, = 2^/2Vpd7,, and Aq ^ 
4.5 eV is the bare gap between the 2p and 3d energy 
levels, unrenormalized by the covalent overlap. 

The elimination of the oxygens in a • • • -0-Ti-O-Ti-- • • 
stripe along the ^/-direction is equivalent to the condition 
VpdcrSy = VpdnSy = 0. Cousequeutly, the absence of the 
local covalency lowers the local symmetry and leads to a 
splitting of the Cg states with the lower energy state cor- 
responding to a 3dsy2_^2 orbital energy with lobes along 
the bond direction and the one at an almost unchanged 
energy of 3dx2_^2 i.e. with lobes in a plane perpendicular 
to the Ti-O^ direction (see scheme of Ti-O^-Ti cluster in 
Fig. [TQ|) with the splitting energies 



pda 



Vp'dJ^ 



A^ 



Ao 



3.2 eV, (3) 

^ -2 eV (4) 



These estimates are for the stripe configuration which 
two O- vacancy neighbors for each Ti. The negative split- 
ting A^^ implies a strong negative shift of the e^-orbital 
with lobes pointing to the vacancy. The e^-orbital is 
shifted below the t2g energy levels, the corresponding 
orbital inversion is observed in the orbital density of 
states in Fig. Efb). In the dimerized Ti-O^-Ti config- 
uration (a) shown in Fig. [TOl the smaller splitting energy 
A^^ ^ — 1 eV results from the missing of only half of 
oxygens in the vacancy stripe configuration. 

The vacancy-induced orbital reconstruction is not re- 
stricted to the LAO/STO interfaces, but has a generic 
character for the titanate surfaces and interfaces. It is 
related to the specific electronic structure of Ti, with 
the covalence-induced splitting between the high-energy 
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(0,0) 



FIG. 11. Band structure for a one uc thick SrTi02 with an 
oxygen vacancy of type (a) in the surface Ti02 layer. The 
thick hues indicate the bands with (a), (b) 3(1^2 _y2 character 
and (c), (d) 3(^3^2 character in the rotated coordinate sys- 
tem shown in Fig. [T]^a). Panels (a), (c) and panels (b), (d) 
display the up-spin and down-spin energies, respectively. 



Cg and low-energy t2g states. For example, in a Ti02- 
terminated (OOl)-surface layer of a SrTiOa slab, the ab- 
sence of the top covalent-bonded oxygens results in the 
five-fold coordination of surface Ti atoms with conse- 
quent d-orbital splitting and reduction of the energy gap 
by about 0.6 eV^^. As a consequence, a surface O- 
vacancy in a Ti02-terminated (OOl)-layer is equivalent 
to the additional condition VpdaSz = in Eq. ([1]), which 
leads to a negative shift of the energy of one of the surface 
e„-states 



The structural relaxation of the vacancy surroundings 
leads to alternating rotations of the z-directed TiiO and 
Ti20 bonds by A9 = ±8° and to a tilting of the reduced 
TiOs-octahedra towards the vacancy center in the {x^z)- 
and (?/,z)-planes (Fig. [TT]) . which increases the overlap. 
Due to the large J > 0, the two-electron local state on 
a Ti-O^-Ti forms a triplet with a triplet energy E2 — J 
where E2 is the two-electron energy in the Ti-O^-Ti clus- 
ter. Due to the planar electron transfer between the 
nearest Ti204-plaquettes, the magnetic triplet state orig- 
inating from the localized 0-vacancies, is spread within 
the interface Ti02-plane and is to be considered as two- 
dimensional magnetic ordering stabilized by the exchange 
splitting of the surface 3d bands. 

The notion that the magnetic ordering at the interface 
has a two-dimensional character is confirmed by calcu- 
lations of the supercells in which the oxygen vacancy is 
residing in more distant layers from the interface Ti02- 
layer. In spite of the occurrence of the local orbital recon- 
struction in the vicinity of the 0-vacancy (Fig. [13]), the 
local magnetic splitting of 3d states is weak and the mag- 
netic moments of Ti in the layer with the vacancies are 
strongly reduced to values of 0.06-0.12 /i^ in the layer sec- 
ond from the interface Ti02-layer, and to 0.006 //b in the 
layer fourth from the interface Ti02. The location of the 
0-vacancies in the layers distant from the interface leads 
to an interlayer exchange between the vacancy-containing 
and interface Ti02 layers, which induces a weak mag- 
netic moment in the interface layer. This effect is similar 
to the double-exchange induced ferromagnetism between 
transition metal ions in different oxidation state^^. Va- 
cancies in the Ti02 layer second from the interface reduce 
the magnitude of the local Ti magnetic moments in the 
vacancy- free interface layer to 0.16 /i^, as compared to 
the magnetic moments in the range of 0.34-0.5 /i^ in 



Ao 



-0.13 eV, 



(5) 



similarly to A^^ of Eq. (j4]) for an interface layer. Our cal- 
culations of a one-uc thick STO-slab indeed confirm the 
formation of a surface 3(i3^2_^2 state in the gap shifted 
by 0.6 eV below the 3dxz band (see Fig. [TT|) . 

In the bulk stoichiometric SrTiOs, the nonmagnetic 
3(i-electron state is associated with a small exchange en- 
ergy splitting of about 0.06 eV separating singlet and 
triplet stated. As follows from the GGA+U electron- 
density contours (Figs. HJ^b) and [6]), the oxygen vacancy 
induces a strong spatial shift of the 3d electron density 
between Ti atoms of about 0.5 A towards the center of 
Ti-O^-Ti-dimer, which increases the overlap exchange in- 
tegrals by about 0.36 eV at the LAO/STO-interfaces of 
type (b) (Fig. [12]). For the dimerized vacancy config- 
uration (a), we obtain from the comparison of the total 
energy of ferromagnetic and antiferromagnetic configura- 
tions the value of the exchange of J = 0.28 eV, however 
we note that the exchange energy is typically overesti- 
mated in the DFT calculations. 
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FIG. 12. Charge density plots across the (OOl)-interface for 
the 0-vacancies of dimer-type (a) and stripe-type (b) in the 
interface Ti02 layer. 
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heterostructures with interface vacancies. 
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FIG. 13. Spin-polarized orbital-projected densities of states 
for the Ti near the O-vacancy(0.5a,0.5a) in the vacancy- 
stripes configuration (b) in the layer second from the in- 
terface Ti02 plane. For this calculation, a supercell 
LAO(4uc)/STO(4uc) has been considered. 



IV. OXYGEN VACANCIES AT THE LAO 
SURFACE 

Another possible source of conducting interfacial 
charge are electrons generated by oxygen vacancies in 
the LAO. Such vacancies are expected to occur predomi- 
nantly in the top AIO2 surface layer— 1^. To explore the 
electronic state of the oxygen vacancies in the AIO2 sur- 
face, we consider first a vacancy placed in the center of 
the (2x1) Al204-plaquette with cy = 1/4 . In this case, 
the surface atomic configuration is described as AI2O3, 
and the charged LAO/STO supercell is expected to be 
doped by the two excess electrons to preserve the overall 
electrostatic neutrality. To study the supercell charging, 
we have determined the two-dimensional charge densities 
in the energy window {Ef — 1 eV; Ep) across the inter- 
face which corresponds to the charge occupation of the 
3(i-conducting bands. 

Fig. [TH presents the charge profile obtained by planar 
integration of the calculated electron densities. The in- 
tegration yields two electron charges per interface unit 
ceh, with 1.25 of electron charges distributed in the sur- 
face and the remaining 0.75 in the STO/LaO-interface 
layer. We note the absence of polar charge produced 
by polar discontinuities in the vacancy- free LAO/STO. 
Consequently, the surface oxygen vacancies suppress the 
polar field of LAO due to the supply of the polarity- 
compensating excess electrons between the surface and 
interface layers. As a result, the interface charging in 
LAO/STO with surface vacancies has a self-doping char- 
acter. 



Fig. [15] shows the projected densities of states for the 
surface Al and O atoms and for the interface Ti. At 
the AIO2 surface, the local states of a mixed A\3sp-02p 
character form in the window {Ep — 1 eV; Ep + 1 eV), 
with the Fermi level located below the top of the local 
band. This implies that the vacancy-doped electrons are 
mobile. In the interface Ti02 layer, the vacancy-doped 
charge occupies the bottom part of the Ti 3dxy conduct- 
ing band, with the distinct ferromagnetic half-metallic 
character reported in Ref. [lo|. A spin polarization is also 
observed in the AIO2 surface states, the surface magnetic 
moments being antiparallel to the magnetic moments of 
the interface Ti02. 

Fig. [T6l shows the spatial distribution of the conducting 
charge, both for the charge which arises from the Al 3sp 
states and for the charge in Ti 3dxy states. The spatial 
distribution of the mobile charges at the AIO2 surface is 
presented in Fig. fTCT a). The conducting charge in the 
Al 3sp-orbitals and O 2p orbitals is shifted towards the 
oxygen vacancies which form stripes along the ^/-direction 
(for the configuration of Fig.[T](b)). In the interface Ti02 
layer, the conducting electrons reside mostly in the 3dxy 
orbitals of Ti atoms (Fig. [TCT b)) which are hybridized 
with 2p states of oxygen atoms. 

For a surface vacancy concentration of cy = 1/8, the 
compensation of the interface polarity by the vacancy- 
released electrons leads to a transition of the AIO2 sur- 
face to the insulating state as can be seen in Fig. [T3 
The insulating surface state occurs due to the transfer of 
the total vacancy-released charge to the interface Ti02 
layer. In the interface Ti02, the transferred charges oc- 
cupy primarily the 3d states of Ti atoms nearest to the 
vacancy. The vacancy-generated electron states at the 
interface are highly spin polarized with Ti magnetic mo- 
ment Mri = 0.56 /iB, enhanced due to a mixed eg-t2g 
character of spin polarization indicated in the orbital- 
projected DOS in Fig. fTTTb). The electron charge profile 
for this case is shown in Fig.[8jb), which demonstrates the 
interface character of the electrons. They are confined 
mostly to the Ti02 layer near the LAO/STO interface. 
The integration of the charge profile in Fig. [81(b) reveals a 
total of two electrons per oxygen vacancy, which implies 
a complete suppression of the interface polarity by the 
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FIG. 14. Charge density profiles along [001] direction for a 
O- vacancy of type (b) located in the AIO2 surface layer in 
the supercell containing 4 uc thick LaAlOs layers and a 4 uc 
thick SrTiOa layer. 
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FIG. 15. Projected DOS (a) for the surface AIO2 and (b) 
interface Ti02 layer in LaA103(4 uc)/SrTi03(4 uc) with one 
0-vacancy of type (b) being present in the AIO2 surface. 



surface oxygen vacancies, a fact supported also by recent 
calculations reported in Ref. 33. Fig. [18] demonstrates 
the calculated x?/-averaged electrostatic potential along 
the [001] -direction. In contrast to the biased macroscopic 
potential in stoichiometric structure (Fig. [9]), the macro- 
scopically averaged potential in the structure with the 
surface vacancies is flat in the L AO-layer, which indeed 
supports the complete suppression of the LAO-polarity 
due to the surface oxygen vacancies. 



V. CONCLUSIONS 

Within DFT-calculations, we considered oxygen va- 
cancies in LAO/STO heterostructures and performed 
studies of the orbital states at the LAO/STO interface al- 
lowing for several types of vacancy arrangements. Using 
generalized gradient approximation (LSDA) with intra- 
atomic Coulomb repulsion (GGA+U), we have shown 
that the oxygen vacancies at the titanate interfaces pro- 
duce a complex multiorbital reconstruction which in- 
volves a change of the occupancy of the Cg states rather 
than of t2g orbitals of the Ti atoms nearest to the oxygen 
vacancies. 

The orbital reconstruction is accompanied by a mag- 
netic splitting of the local Cg and interface dxy orbitals. 
This reconstruction generates a two-dimensional mag- 




FIG. 16. Charge density plots for (a) the AIO2 surface and (b) 
interface Ti02 layers with the O- vacancy in the AIO2 surface. 
Here a supercell LAO(4uc)/STO(4uc) has been considered. 



netic state not observed in bulk SrTiOa. Moreover, oxy- 
gen vacancies placed in the Ti02 layer farther away from 
the interface induce a sizable magnetic moment only in 
the Ti02 interface layer. Also, oxygen vacancies in the 
AIO2 surface where they are expected to be in the most 
stable configuration, generate a magnetic moment only in 
the interface titanate layer. In this latter case, the elec- 
tronic reconstruction mechanism due to polar catastro- 
phe is suppressed for a vacancy concentration of cy > 1/8 
by the charge introduced by vacancies. The surface is 
then free of charge carriers which implies the formation 
of an insulating state at the AIO2 surface. 



In configurations with vacancy stripes, we have found 
an orbital separation of the charge and spin degrees of 
freedom with the vacancy-released charge carriers local- 
ized in e^-orbitals, and the spin polarization occurring 
predominantly in t2^-states. Moreover, we have provided 
evidence for the generic character of the two-dimensional 
magnetic state at titanate surfaces and interfaces. 
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FIG. 17. Projected DOS (a) for the surface AIO2 and (b) 
interface Ti02 layers in LaA103(4u.c.)/SrTi03(lu.c.) with 
one 0-vacancy per eight 0-atoms in the AIO2 surface. The 
positions of Ti, Al, La and O atoms have been fully relaxed 
in the (x, 2/)-planes and in the ^-direction. 
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FIG. 18. Profile of x2/-averaged electrostatic potential along 
the [001] direction with cv = 1/8-concentration of surface 
oxygen vacancies. The black and red profiles correspond to 
the x?/-averaged and macroscopically averaged potentials, re- 
spectively. 
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